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We study spinless fermions in a flux threaded AB2 chain taking into account nearest-neighbor 
Coulomb interactions. The exact diagonahzation of the spinless AB2 chain is presented in the limit- 
ing cases of infinite or zero nearest-neighbor Coulomb repulsion for any filling. Without interactions, 
the AB2 chain has a fiat band even in the presence of magnetic flux. We show that the respective 
localized states can be written in the most compact form as standing waves in one or two consecutive 
plaquettes. We show that this result is easily generalized to other frustrated lattices such as the 
Lieb lattice. A restricted Hartree-Fock study of the V/t versus filling phase diagram of the AB2 
chain has also been carried out. The validity of the mean-field approach is discussed taking into 
account the exact results in the case of infinite repulsion. The ground-state energy as a function of 
filling and interaction V is determined using the mean-field approach and exactly for infinite or zero 
V. In the strong-coupling limit, two kinds of localized states occur: one-particle localized states 
due to geometry and two-particle localized states due to interaction and geometry. These localized 
fermions create open boundary regions for itinerant carriers. At filling p = 2/9 and in order to avoid 
the existence of itinerant fermions with positive kinetic energy, phase separation occurs between a 
high-density phase (p = 2/3) and a low-density phase (p — 2/9) leading to a metal-insulator tran- 
sition. The ground-state energy reflects such phase separation by becoming linear on flUing above 
2/9. We argue that for filling near or larger than 2/9, the spectrum of the t-V AB2 chain can be 
viewed as a mix of the spectra of Luttinger liquids (LL) with different fillings, boundary conditions, 
and LL velocities. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.10.Pm, 75.10. Pq 



I. INTRODUCTION 

The field of itinerant geometrically frustrated elec- 
tronic systems has attracted considerable interest in the 
last two decadeJilt^Il A simple example of a geometri- 
cally frustrated lattice is the AB2 chain, also designated 
by diamond chain or bipartite lozenge chairPl^. This 
is a quasi-one-dimensional system consisting of an one- 
dimensional array of quantum rings. Such system can be 
used to model ML2 (metal-ligand) chainP2] and azurit^ 
and can be generalized to molecular systems displaying 
similar topologj^. Also, given nowadays nanofabrica- 
tion techniques, such as electron beam lithography, a di- 
amond chain-like system can, in principle, be built from 
scratclP. A particular feature of the band structure of 
these frustrated systems is the presence of one or several 
flat bands in the one-particle energy dispersion, which 
reflect the existence of localized eigenstates of the geo- 
metrically frustrated tight-binding HamiltoniarPtl^. 

The ground states of the Hubbard model for an AB2 
geometry are well studiecPMSl as well as for other frus- 
trated lattices such as the sawtooth chairf^I, the Kagome 
chaiiPi^, etc. Some of these frustrated lattices fall onto 
the category of the cell construction lattices proposed 
by TasakPl or the category of line graphs lattices pro- 
posed by Mielk^. The approach followed in some of 
these studies relies in the fact that the Hubbard interac- 
tion is a positive semidefinite operator and is limited to 
cases where the lowest band is a fiat band or the chemi- 



cal potential is fixed at the flat-band energ jff. The frus- 
trated systems have usually been studied in the absence 
of flux'^ * In the case of AB2 Hubbard chain, the flux 
dependence of the ground-state energy has been studied 
but again only for chemical potential fixed at the fiat 
band energy^. Many different ground states are possi- 
ble in the AB2 Hubbard chain, leading to a great variety 
of properties such as fiat-band ferromagnetism or half- 
metallic conduction, depending on the values of filling, 
interaction or magnetic fielcP'. 

The case of spinless fermions in a AB2 lattice, tak- 
ing into account nearest-neighbor Coulomb interaction, 
is simpler than the Hubbard model due to the absence 
of spin degrees of freedom. We will designate the respec- 
tive Hamiltonian for such system by t-V AB2 Hamilto- 
nian. The t-V model, in its strictly one-dimensional ver- 
sion (ID), can be mapped into the anisotropic Heisen- 
berg model (more precisely, the XXZ or Heisenberg-Ising 
model) by the Jordan- Wigner transformatiorPSl, whose 
Bethe ansatz solution has long been knowrP^. For quasi- 
ID models such as the one discussed in this manuscript, a 
Bethe ansatz solution is not possible. However, a Jordan- 
Wigner transformation into the XXZ AB2 chain should 
be possible using the extension of the Jordan- Wigner 
transformation to two dimensions which has been dis- 
cussed by several authors^^^^^. In the case of a square 
lattice, this transformation requires the introduction of 
a gauge field which in contrast to the one-dimensional 
case affects the energy spectruirl^Si. In particular, the 
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strong-coupling limit of the repulsive t-V model un- 
der a Jordan- Wigner transformation is mapped into the 
strongly anisotropic antiferromagnetic Heisenberg model. 
During the last years, many studies of the antiferromag- 
netic Heisenberg model in geometrically frustr ated lat- 
tices have been carried ouP^IMnillllMllIlIlini Under 
high magnetic fields but below the saturation field, the 
ground states of these frustrated magnetic systems con- 
sist of localized and independent magnons created in a 
ferromagnetic background. As a consequence of these lo- 
calized magnons, quantized plateaus have been found in 
the respective magnetization curve^^^. 

In this manuscript, the exact diagonalization of the 
spinless AB2 chain is presented in the limiting cases of 
infinite or zero nearest-neighbor Coulomb interaction for 
any filling and in the presence of magnetic flux. Without 
interactions, the AB2 chain has a flat band even in the 
presence of magnetic flux. A simple construction of the 
localized states that generate the flat bands both in the 
presence and absence of flux is presented and generalized 
to arrays of AB„ quantum rings. The flat band gener- 
ates a plateau in the ground-state energy as a function 
of filling (for fillings between 1/3 and 2/3). A restricted 
Hartree-Fock study of the V/t versus filling phase dia- 
gram has been carried out. For finite V , the mean-field 
ground-state energy increases in relation to the indepen- 
dent fermions ground-state energy, but remains negative 
for fillings lower than 2/3. For filling larger than 2/3, it 
becomes approximately linear reflecting the existence of 
nearest-neighbor occupied sites. In the mean-field phase 
diagram, a uniform density phase is found at low filling. 
For filling larger than 1/3, one of the mean-field solutions 
for the density difference between A sites and B and C 
pairs of sites disappears (due to Pauli's exclusion princi- 
ple) . The validity of the mean-field approach is discussed 
taking into account the exact results in the limiting cases. 

The strong-coupling limit is particularly interesting 
due to the presence of two kinds of localized states: one- 
particle localized states due to geometry and two-particle 
localized states due to interaction and geometry. Local- 
ized states may lead to flat bands or mix with itiner- 
ant states creating open-boundary regions for itinerant 
carriers. The localized fermions due to geometry that 
appear in the frustrated AB2 chain have a direct corre- 
spondence with the independent localized magnons in the 
near-saturation magnetically frustrated systems. Note, 
however, that the equivalent to the two-particle localized 
fermions in the frustrated AB2 chain has never been men- 
tioned in studies of the near-saturation antiferromagnetic 
Heisenberg model in geometrically frustrated lattices, as 
far as we know. The zero-temperature phase diagram 
obtained plotting the ground-state energy as a function 
of fllling displays an interesting quantum critical point 
at filling p = 2/9 where a metal-insulator transition oc- 
curs. This transition refiects a phase separation between 
a high density phase [p = 2/3) and a low-density phase 
(p — 2/9) that occurs at fillings larger than p — 2/9. 
It is worthwhile to emphasize that phase separation is 



know to occur in the attractive t-V chairpiE3 but it does 
not occur for repulsive nearest-neighbor interactions. In 
the t-V AB2 chain, we have shown that phase separation 
occurs for strong repulsive nearest-neighbor interactions. 
The particularity about the filling p = 2/9 is that the 
reduced effective lattice (infinite nearest-neighbor repul- 
sion effectively reduces the number of available sites for 
fermions) becomes half-filled at p — 2/9. Phase separa- 
tion occurs in order to avoid the existence of itinerant 
fermions with positive kinetic energy. We also discuss 
whether the strong coupling t-V AB2 chain is a Luttinger 
liquid. We argue that while for low filling, the low energy 
properties of t-V AB2 chain can be described by the spin- 
less Luttinger Hamiltonian, for filling near or larger than 
2/9, the AB2 set of eigenstates and eigenvalues becomes 
a complex mix of the sets of eigenstates and eigenvalues 
of LLs with different sizes, fillings, boundary conditions, 
and LL velocities. 

The remaining part of this paper is organized in the 
following way. In Sec. II, the model is defined. In Sec. 

III, the eigenstates of the AB2 tight-binding model in 
the presence of magnetic flux are found and the ground- 
state energy as a function of filling is determined for any 
flux. A generalization of the flat-band states to AB„ 
chains or even more complex lattices is discussed in Sec. 

IV. In Sec. V, the study of the restricted Hartree-Fock 
phase diagram of the AB2 chain for V/t versus fllling 
is presented. In Sec. VI, the strong-coupling limit of 
the model is addressed. The Luttinger liquid description 
of the strong coupling t-V AB2 chain is also discussed. 
In Sec. VII, we show that the results obtained in the 
previous sections are relevant for the extended Hubbard 
model in the same geometry. In Sec. VIII, we conclude. 



II. THE AB2 CHAIN 

In Fig. [ija) , a diamond ring is shown with a magnetic 
flux (j) threading each diamond plaquette and a mag- 
netic fiux (j)i threading the inner ring. The inner sites 
in Fig. [TJa) are denoted as C sites and the outer sites as 
B sites. Note that there are two ways to close the linear 
AB2 chain shown in Fig. [ijb), either by leaving the B 
sites or the C sites in the interior of the ring. The two 
situations are physically equivalent and therefore we will 
assume that we have closed our ring so as to leave the 
C sites as inner sites. The system may be pictured as 
two rings, an outer one and an inner one, as shown in 
Fig. [l](a), so that an electron traveling through the outer 
ring sees an effective fiux (po while and electron traveling 
through the inner ring sees a different effective flux 4n. 
It will prove itself useful to introduce an auxiliary flux 0' 
such that 

0o = <^' + Afct 



FIG. 2. The Hamiltonian of the diamond chain threaded by 
an arbitrary flux can be mapped into a non-interacting peri- 
odic Anderson-hke model with a basis of three sites per unit 
cell and with hopping amplitudes controlled by magnetic flux. 
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FIG. 1. (a) The AB2 chain consists of an one-dimensional 
array of diamond rings. Under magnetic field, each diamond 
ring is threaded by a magnetic flux <j) and the inner star is 
threaded by a flux (f>i. This system can be pictured as an 
outer ring with flux (j)o and an inner ring with a flux cj>i. (b) 
The diamond chain when </)' = 0. 



where is the number of cells of the diamond ring. In 
the case of zero both effective fluxes still remain non- 
zero if <j) is nonzero. In this case, the Peierls phase can be 
equally distributed in each plaquette, restoring the trans- 
lational invariance in each ring as shown in Fig. [TJ^b) . For 
= or TT, the lattice is invariant in the "flip" of one pla- 
quette (so that B and C sites exchange places) , reflecting 
a local Z2 symmetrjl2ll. 

Considering nearest-neighbor Coulomb interactions 
and particle hoppings, the t-V Hamiltonian for an AB2 
chain with N(, unit cells (or diamonds) is 



H = Ho + V 
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(2) 



where V is the value of the interaction and 

+ e-'^'/^^=(C]A,+A]+iQ) 



(3) 



H.c. 



Here we have chosen a gauge such that the Peierls phases 
are equally distributed in the inner ring and in the outer 
ring. 



III. TIGHT-BINDING LIMIT 

For simplicity, we will assume a flux configuration such 
that </)' = in this section. The general case is considered 
at the end of this section. The tight-binding Hamiltonian 
for AB2 chain with Nc cells is 



j 

+ A]^, {e-'^/^Bj + e'-^I^C,) + H.c. 



(4) 



Employing the transformations 




g-i0/4 gi0/4 




(5) 



the system may be mapped into an Anderson-like model 
as depicted in Fig. [2] Its Hamiltonian is given by 
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H.c. 



(6) 



where 



h = V2t, (7) 
t2 = V2cos{(f)/2)t, 
h = v/2sin(0/2)t. 

We note that the possibility of mapping a chain Hamilto- 
nian into a periodic Anderson model has been already ex- 
plored by Gulacsi et al by mapping the triangular-chain 
Hamiltonian into a periodic Anderson model'^'*. In our 
case, when the flux is either zero or tt, one has localized 
eigenstates, which lead to flat bands in the energy dis- 
persion (see Fig. [4| . For zero flux the system becomes 
a tight-binding ring of A and b sites with independent 
c sites (see Fig. [2| , having plane- wave eigenstates in the 
ring and localized eigenstates at c sites, as depicted in 
Fig. [sjja). On the other hand, for = tt, the system 
behaves as a set of Nc independent systems with three 
sites and all the eigenstates are localized as depicted in 
Fig. [sf^bj^. The eigenvalues for an arbitrary value of 
flux are given by 




(a) 




(b) 

FIG. 3. Eigenstates of the tight-binding AB2 chain for (a) 
<j!> = and (b) (f) — n. For zero flux, itinerant (red dashed 
line) as well as localized eigenstates (blue dotted line) are 
present. For (j) = tt, all states are localizecP^. 
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FIG. 4. Dispersion relation of the diamond chain for different 
flux values (we have omitted the flat band = since it 
is not affected by the flux). For (p = n, all bands are flat 
reflecting the local Z2 symmetry of the Hamiltonian. 



Eflat = 0, 
e± = ±2ty/l + cos{(l)/2) cos(fc). 



(8) 



In Fig. |4j the dispersion relation for several values of flux 
is shown. The flat band e/c = is omitted since it is not 
modifled by the presence of flux. We introduce here the 
fermionic operators fl , aj, and dj^ , which create electrons 
in the single-electron bands, more precisely, fl creates a 
particle with momentum k on the top band, aj, on the 
flat band, and on the bottom band. The eigenstates 
creating operators are given by linear combinations of the 
creation operators on sites A, B and C, and for the flat 
band, the expression of the creation operator is rather 
simple. 



4 = 



y/l + COS (0/2) COs(fc) 



COS (</./4 + k/2) Cl 



cos ((/./4 - k/2) bI 



(9) 



As can be concluded from the above results, the dia- 
mond chain presents always a flux independent disper- 
sionless band and is gapless for zero flux. At flnite flux a 
gap opens between the bo ttom and the top bands. This 
gap is given by Ae = 4t-\/l — cos(0/2), while the band- 
width of the system is = 4^^^! -I- cos {4>/2). 

InFig.|4) it can be observed that for k — tt/2, two flux 
independent eigenvalues appear. This condition is only 
fulfilled if the number of unit cells is a multiple of four, 
Nc — An with n integer. We should also note that the 
system is 27r periodic on the flux when N^. is even and 47r 
periodic when Nc is odd. A flux of 27r interchanges the 
energies of the dispersive states with k = Q and k = tt. 
These two values of momenta are allowed when Nc is 
even. The value fc = tt is forbidden when Nc is even and 
the period of the system is An on the flux. 

The eigenvalues for an arbitrary value of flux 0' are ob- 
tained following similar steps and the energy-dispersion 



relations become 

fflat = 0, 



= ±2^1 + cos((/)/2) cos{(t)'lNc + k). 



(10) 



The insertion of an extra flux (j)' in the AB2 chain only 
translates the energy-dispersion relation by (j)' l^c and 
therefore does nothing to the flat band. This result 
is expected since it is similar to what is observed in 
ID quantum rings. One can use (j)' to control the mo- 
mentum for which the top/bottom band reach its maxi- 
mum/minimum energy. 



A. Density of states, filling and ground-state 
energy 

Let us consider the thermodynamic limit (fc contin- 
uous) and calculate the density of states (DOS) of the 
tight-binding AB2 chain. The DOS of the flat band is a 
Dirac-delta function at e = 0. The dispersive bands are 
one-dimensional bands and the combined DOS of the two 
dispersive branches can be written as 



1 

Vt 



y^cos2(0/2) - (e'2 - if 



(11) 



when e belongs to the energy intervals associated with 
the dispersive bands and where e' — e/2t. The full DOS 
of the system is therefore, D{e) — D±{e) -f 5{e). 

Filling is deflned as the number of electrons N per site, 
p — N/Ng — N/3Nc, where N^ is the number of sites. 
Due to the symmetric nature of the energy spectrum and 
since each band can accommodate one electron per unit 
cell, we know that half flUing occurs for Ep — 0. Since 
the flat band contains Nc states then when p E [1/3, 2/3] 
one has Ep = 0. If Ep lies on the bottom band we have 
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P = 



D_(e)d6, 



(12) 
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FIG. 5. Ground state energy in the thermodynamic hmit as 
a function of filhng for several values of flux. Eqs remains 
constant when filling the flat band, p G [1/3,2/3] achieving 
its minimum value in that interval for (j> — tt. For = tt, _Egs 
has linear behavior in p, since all bands are flat. 



where eb^mm — — 2\/2i cos(0/4) is the bottom of the band 
and where the factor of 1/3 is due to the fact that the 
integral of the DOS over a band is one and each band 
contributes equally to the DOS. Using Eq. (11 1 we have, 



P = 



1 



arctan 



b,min 



arctan 



Ycos2(^/2P^ 

ycos^Wsy^^y 



(13) 



where 



2t ) 



1, 



Ef 
2t 



(14) 



- 1. 



If, on the other hand, Ep lies on the top band, one can use 
the symmetry of the density of states to write p = 1 — p/i, 
where 



1 



Ph = 



- Ep 



^-(e)de. 



(15) 



is given by Eq. (13). 

Given the preceding results, we are able to calculate 
the ground-state energy of the system as a function of 
filling. The ground-state energy is given by 



i:>(e)ede. 



(16) 



while Eqs. (13) and (14) give us the relation between 



the Fermi energy and the filling. For a general value of (j) 
there is no simple analytical expression and the ground- 
state energy should be calculated numerically. However, 
when = TT it is very simple to derive an exact result. 



given that the system only has three energy levels, e = 
0, ±2t and that each can take Nc electrons. In this case, 
the ground-state energy per unit cell while the bottom 
band is not fully filled is given by E'qs = —2tN/Nc = 
—dtp. At p = 1/3 the bottom band gets fully filled and 
from p = 1/3 to p = 2/3 we are filling the middle band 
whose energy is e = 0. Therefore, Eqs = — 2t, p e 
[1/3,2/3]. For p > 2/3 we start filling the upper band, 
whose energy is e = 2t. In this situation the ground-state 
energy is given by Eqs = 6t(p — 2/3) — 2t. 

In Fig. [5]we plot the ground-state energy as a function 
of filling for several values of flux. As can be seen for 
p e [1/3,2/3] the ground-state energy remains constant 
since we are filling the flat band = 0. The ground- 
state energy is obviously even around = 0.5 due to 
the symmetric nature of the dispersion relation. We also 
see that for p ^ [1/3, 2/3] the dependence of the ground- 
state energy on the filling departs from its linear behavior 
when (f) ^ TT. 



IV. FLAT BANDS AND LOCALIZED STATES 
FOR ARRAYS OF QUANTUM RINGS 

It is well known that the single-particle flat band eigen- 
states of a geometrically frustrated lattice can be written 
as a set of localized eigenstates which are translated ver- 
sions of the same state | V'loc) ■ This single-particle state is 
non-zero in a small lattice region and it is an eigenstate of 
the tight-binding Hamiltonian (a zero energy eigenstate 
in the case of the AB2 chain). This means that if we write 

the eigenstate as \tp\oc) = J2i '^i'^ll^) ' where i runs over 
all lattice sites such that 7^ 0, then t^^ai = ear = 0, 
where r can be any site in the lattice such that = 0, 
tri is the hopping constant between sites r and i and e is 
the kinetic energy of the state. In the particular case of 
AB2 chains, e = and the restrictions on sites i and r 
can be lifted, i.e., J^i^ri^i = ear = 0, where i runs over 
all lattice sites and r is any site of the lattice. We pro- 
pose a particular perspective (but equivalent) for the case 
of any array of quantum rings similar to the AB2 chain 
which easily allows us to identify the localized states in 
the absence or presence of flux. Let us also assume for 
now that the array is such that the shared sites between 
consecutive rings are directly opposite, and therefore the 
number of sites in each ring is even. 

For zero flux, (j> = (po = ipi = 0, all energy levels of 
the tight-binding ring are doubly degenerate with ex- 
ception of the lowest and highest energy levels (A; = 
and k = tt) and the respective eigenstates have opposite 
momenta, reflecting the time-reversal symmetry of the 
Hamiltonian. In the particular case of the AB2 chain, 
the respective ring has four sites and there is only one 
degenerate level corresponding to zero energy. Subtract- 
ing the states of opposite momenta, one obtains a stand- 
ing wave with nodes at i = and i — Nr/2 (ring sites 
are numbered clockwise from i = to i = Nr — 1). For 
Nr > 2n with n > 2 one has more than one flat band. 
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More precisely, one has N^-mg/'^ — 1 flat bands with en- 
ergies e„ = —2t cos{2Trn/ Nr) with n — 1, . . . , Nr/2 — 1. 
Since these nodes coincide with the sites shared by con- 
secutive rings, one can say that for these standing wave 
states, the ring becomes uncoupled to the rest of the 
chain. Furthermore, if Nr > 2n with n > 2 the standing 
waves may have more than two nodes and more complex 
lattices can have such localized states. Note however that 
the nodes position must be an integer. We conclude that 
in order to find flat band eigenstates of arrays of quan- 
tum rings, one may construct standing waves such that 
the nodes coincide with sites shared between different 
quantum rings. For instance, for Nr = 8, a square lattice 
is possible as shown in Fig. |8] 

These arguments agree with the results obtained for 
the AB2 chain. For (f> = (j),, = 4>i = 0, one has a\. = 
[bI-CI]/V2, which can be Fourier transformed to obtain 
al = [Bj - C\\l\f2. The single-particle state a|'|0) is a 
localized eigenstate in cell i with zero energy. 

Let us consider now = but 0o = ~4>i 7^ 0. 
Then the tight-binding constants in one ring become 
j+i = e-^-^o/^^-'-Ji, for j = 0, . . . , 1 and t^j+i = 

e'''^°/2We-Ji^ for i = Nr/2, ...,Nr - I so that the total 
Peierls phase is zero reflecting the fact that i/i = 0. A 
simple gauge transformation 

c] ^ e'^-^c], J=0,...,^ 

c] ^ e'^-(^'-^")c], j = i^ + l,...,Ar^-l 

restores the translation invariance with a zero Peierls 
phase and the previous construction of localized states 
can be applied. Therefore, the same flat bands will be 
present. In the case of the AB2 chain, the flat-band states 
will have the form aj|0) = e-*5fe[E| - C\]/^/2\Q), but 
the phase term is obviously irrelevant. 

Let us address now the case of </> ^ but </)' = (the 
case with 7^ is obtained following the same steps after 
a gauge transformation as explained for the = case) . 
When 7^ 0, eigenstates of the tight-binding ring with 
opposite momenta k and — fc do not have the same energy 
(see Fig. [6]) and therefore the standing wave obtained 
from the subtraction of these states is not an eigenstate 
of the tight-binding Hamiltonian of the ring. One can 
overcome this difficulty considering two consecutive rings 
and following a certain path in those rings such that the 
total Peierls phase in that path is zero as shown in Fig. [7] 
If one constructs standing waves for a ring of 2Nr sites 
with the additional condition of extra nodes at sites Nr/2 
and 2Nr — Nr/2 (sites A and B in Fig. [t]), this state can 
be folded (making site A overlap with site B) to give an 
eigenstate of the two rings threaded by flux with nodes 
at the extremities. Note that the larger ring has non- 
zero Peierls phases associated with each hopping but the 
total phase is zero. So after a gauge transformation to 
eliminate these phases, the states of opposite momenta 
are degenerate and the standing wave is an eigenstate of 
the Hamiltonian of the two rings. Again, one or more flat 
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FIG. 6. Dispersion relation of a tight-binding ring with eight 
sites: (a) without flux; (b) threaded by a magnetic flux (p = 
7r/5. For zero flux, states k and ~k are doubly degenerate 
(with the obvious exception of fc = and fc = tt). When 
7^ 0, the eigenstates of the tight-binding ring with opposite 
momenta fc and — fc do not have the same energy reflecting 
the loss of time-reversal symmetry. 
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FIG. 7. In the presence of flux, considering two consecutive 
rings and following flrst the continuous path (in red) and then 
the dotted path (in black) the total Peierls phase in that path 
is zero (the path is clockwise for the left ring and counter- 
clockwise for the right ring). If one constructs standing waves 
for a ring of 2Nr sites with the additional condition of extra 
nodes at sites A and B, this state can be folded to give an 
eigenstate of the two rings threaded by flux with nodes at the 
extremities and at site A. 

bands are possible depending on the size of the quantum 
ring. 

The gauge transformation for the larger ring will be 

c] ^ e^^-(^-^)c], J = ^f + l,...,Nr 

c] e-^*-(2^'--'\], J=Nr + ^f + l,...,2Nr~l 

where the sites have been numbered clockwise in the 
larger ring. In the case of the AB2 chain, after the 
gauge transformation the localized state (not normal- 
ized) will be (c| — C3 + C5 — Cy)|0), which, invert- 
ing the gauge, transformation corresponds to the state 

[e-^^(c| - 4) + e'^^{cl - 4)]|0)- This state writ- 
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FIG. 8. Localized states in the case of the Lieb lattice with 
or without magnetic flux. The introduction of magnetic flux 
extends the localized state to two plaquettes. These plaque- 
ttes are equivalent to the quantum rings discussed in the text 
and the localized state is a standing wave with nodes in the 
plaquettes vertices. 



ten in terms of the operators of the AB2 chain becomes 

[e-*^(B]-C]+i)+e'^(B]+i-C])]|0) or equivalently 

[e-'i [b] - C]_^i) + e*^ [b]^^ - C])] |0) . These are exactly 
the localized eigenstates obtained Fourier transforming 
Eq. (|9|, 



2e^'=/2 ^1 + cos {(j)/2) cos(fc)4 



p-«(</'/4-fc)]^t 



,i(0/4+fc) 



-i0/4l/-(t 



(17) 



Note that these localized states may overlap in real space, 
that is, they constitute a basis of the subspace of localized 
states but not an orthogonal basis. Such impossibility of 
constructing orthogonalized Wannier states for certain 
lattices with flat bands under a uniform magnetic field 
has been mentioned befor^^H 

As an example of application of the previous argu- 
ments, we show in Fig.|8]the localized states for the Lieb 
lattice. It is known that this lattice displays a flat band 
for zero and finite magnetic flux. Aoki et al. have found 
the respective localized state for zero flux as well as by 
inspection a localized "elongated ring state" for finite 
flusP^. The zero fiux localized state agrees with the one 
shown in Fig.jsj but our localized state for finite flux 
shown in Fig. [8j is considerably more compact than the 
"elongated ring state" of Ref.l^. We expect that the pre- 
vious arguments can be applied to other lattices that fall 
onto the Lieb lattice category, that is, bipartite lattices 
with different numbers, ua and ns, of A and B sublattice 
sites in a unit celP^. Another example of simple appli- 
cation is the Be class superhoneycomb lattic^2l! Again 
our argument justifies the fact that the flat band in this 
system remains flat for finite magnetic field. 



V. MEAN-FIELD RESULTS FOR GENERAL V 

In this section, we present a mean-field study of the t-V 
AB2 chain taking into account nearest-neighbor Coulomb 



interaction. The results obtained must be interpreted 
with caution, having in mind the known drawbacks of this 
approach. One of these drawbacks is the fact that this 
approach neglects thermal fluctuations, which as stated 
by the Mermin- Wagner theorem prevent long range or- 
der at finite temperature. In this section, only the zero- 
temperature case is addressed, so one avoids this prob- 
lem. Even at zero temperature, the mean-field approach 
overestimates the existence of an ordered phase, since 
quantum fiuctuations oppose the emergence of an or- 
dered phase and these quantum fiuctuations are particu- 
larly strong in quasi-lD systems. Note however that the 
ground state of quasi-lD systems may in fact be ordered 
despite these quantum fluctuations. 

When considering nearest-neighbor Coulomb interac- 
tions, we will be interested on the density of particles at 
A, B and C sites. We assume the particle density on site 
X to be the average number of particles per number of 
unit cell, px = N x /A^c- We also assume that the particle 
density at B and C sites is the same. Let p denote the 
total particle density. 



PA + 2pi 



(18) 



where pA and pB are the particle densities at A and B 
sites respectively. We then have 



.PB, 



(19) 



In this situation, the interaction part of the mean-field 
Hamiltonian can be written as 

iJint = 2y ^ [2pBnf -f pAuf + pAuf] ~ AN^PAPB- 

(20) 

While there is no simple expression for the mean-field dis- 
persion relation for general 0, a simple expression exists 
for zero flux, 



Efiat = "iVpA, 



e±^V (pa + 2pb)± JSt^ cos2(A:/2) + A^., 



(21) 



where Ay — V (pA — '^Pb)- Again, a flat band is present, 
but its energy level depends on the density at sites A. 

We define the order parameter as the excess of density 
at the A sites. 



Ap: 



PA - "^PB 
PA + ^Pb 



(22) 



Due to the equivalence of the B and C sites, we re- 
gard our system as being a dimerized system consist- 
ing of two alternating types of sites: A sites and BC 
pseudo-sites. Note however that this picture is to be in- 
terpreted with caution since BC pseudo-sites are not real 
sites and can accommodate twice as many electrons as A 
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P 

FIG. 9. V jt versus filling mean-field phase diagram for the 
spinless diamond chain considering nearest-neighbor Coulomb 
interactions and TVc = 128. For p < 1/3 we can have a uniform 
density phase (Ap = 0) or a phase with excess of density at 
A or BC pseudo-sites (Ap > or Ap < 0, respectively). 
For p > 1/3, Pauli's exclusion principle breaks the symmetry 
between A and BC pseudo-sites. A uniform density phase 
is not possible anymore and the density of particles can no 
longer be situated only at A sites. It can however, for p < 2/3, 
be situated only at BC pseudo-sites. For p > 1/3, due to 
Pauli's exclusion principle, the density of particles is required 
to be spread among A and BC sites and for p = 1 the order 
parameter is Ap = —1/3, which implies a uniform density of 
particles between the real A, B and C sites. 



sites. One must note that in a general situation we have 
Ap G [— Ijl]- Since < pa,Pb < 1, then for p > 1/3, 
Ap is hmited to the interval [— 1, Apmax], where Ap^ax 
Ues in ] — 1,0[ and whose value decreases with increas- 
ing p. For the same reason, for p > 2/3, Ap is limited to 
] — 1, —1/3], and not only the upper limit of Ap decreases 
with increasing p, but also the lower limit of Ap increases 
with increasing p. For p = 1 we have Ap = —1/3, cor- 
responding to an equal density of particles on every site 
(the only possible state when the system is completely 
filled). 

The phase diagram of the system is depicted in Fig. 
[9j a uniform density phase can exist only for a filling 
p < 1/3. For p < 1/3, starting in the uniform density 
phase, by increasing the interaction we are able to lo- 
calize an excess of electron density at A sites or at BC 
pseudo-sites, both situations being symmetric. By fur- 
ther increasing the interaction we are able to localize all 
the fermions at A sites or at BC pseudo-sites where both 
situations remain symmetric. As a consequence the order 
parameter Ap can be double valued for these values of 
filling. In the region 1/3 < p < 2/3 one no longer can lo- 
calize the full electrons density at A sites while one can at 
BC sites and therefore although the Hamiltonian treats A 
and BC pseudo-sites equally, Pauli's exclusion breaks the 
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FIG. 10. Order parameter as a function of filling for V/t = 
0.5, 3 and 20 and Nc = 200. The changes of slope of the 
order parameter indicate the phase transitions. Note that the 
absence of a positive solution of Ap for p > 1/3 reflects the 
fact that the A and BC pseudo-sites may be occupied by one 
and two fermions, respectively. As a consequence, for p < 1/3 
the order parameter can be double valued since we are able 
to localize all fermions at A sites (Ap > 0) or at BC pseudo- 
sites (Ap < 0). Since for p > 1/3 we can no longer localize 
all fermions at A sites, the order parameter can no longer be 
positive. 
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FIG. 11. Ground state energy per site as a function of p for 
V/t = 0, 1, 3 and 20 and Nc — 200, obtained in the mean- 
field approach. For finite V, no flat region appears and the 
minimum energy is shifted to lower flUing. The ground-state 
energy remains negative for p < 2/3. For p > 2/3, the ground- 
state energy is almost linear in the flUing with a large slope, 
since nearest-neighbor pairs are being created. 

symmetry between A and BC pseudo-sites, lowering the 
symmetry of the system. Consequently, for p > 1/3, the 
order parameter can no longer be double valued. In the 
region p > 2/3, again due to Pauli's exclusion principle, 
one is not able to fully localize the density of electrons 
even at BC pseudo-sites and only one phase remains. 

In Fig. [TOj the order parameter as a function of filling 
is shown for V/t = 0.5, 3 and 20 and Nc = 200. The 
changes of slope of the order parameter indicate the phase 
transitions. The absence of a positive solution of Ap for 
p > 1/3 reflects the fact that the A and BC pseudo-sites 
may be occupied by one and two fermions, respectively. 
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In Fig. [TTJ we show the ground-state energy per site as 
a function of p for V/t = 0, 1, 3 and 20 and = 200, 
obtained in the mean-field approach. For finite V, no 
flat region appears and the minimum energy is shifted 
to lower filling. Recalling Eq. (21), one concludes that 
a constant term 3Vp^ is added to the non-interacting 
ground-state energy when Ap = 0, which shifts the min- 
imum energy to lower p when V is small (Ap = for 
low p). The ground-state energy remains negative for 
p < 2/3. For p > 2/3 and large V, the ground-state 
energy is almost linear in the filling with a large slope, 
since nearest-neighbor pairs are being created. 



VI. THE ABa CHAIN IN THE 
STRONG-COUPLING LIMIT ^ oo 

Making t/V a small parameter, one drives the t-V AB2 
model into the so-called strong-coupling limit. There are 
two equivalent ways to reach this limit, either increas- 



ing V or reducing t. li t = 0, the fermions are localized 
and all states with the same number of pairs of nearest- 
neighbor occupied sites, + + "■F')' 
degenerate. This degeneracy is much lower compared 
with the ground-state degeneracy of the — > 00 Hub- 
bard AB2 model. This degeneracy is lifted if t/V is finite 
and up to first order in t, the eigenvalues are obtained 
diagonalizing the Hamiltonian within each of the degen- 
erate subspaces. The Hamiltonian within each subspace 
is obtained using the Gutzwiller projection operators Pi, 
which project onto the subspace with I pairs of nearest- 
neighbor occupied sites. The set of eigenstates and eigen- 
values of Z = subspace of this model can be determined 
as we will show below, relying at least for some of the 
eigenstates in the knowledge of the solution of the t-V 
chaiiP3. 

In order to simplify the comparison between the t-V 
chain and the AB2 chain, we number the sites in a differ- 
ent way, using odd numbers for A sites and even numbers 
for B and C sites. 



Let us also define the operator = (1 — c^q) where c can be an operator A, B or C . In the strong-coupling limit 
V/t » 1, we obtain for the ground state subspace {I = 0) the projected Hamiltonian with cj) = Q (but (f)i = (j)o ^ 0) 



j odd 



i e Vi 



i e Vi 



H.c. 



where Vj is the set of sites nearest-neighbors of the pair of sites j and j + 1 (excluding these sites) and the product 
of hole occupation numbers reflects the condition that a nearest-neighbor pair of occupied sites is not created when 
a particle hops. 

I 



Again, we emphasize that states belonging to subspaces 
with I ^ pairs of nearest-neighbor occupied sites will 
be discarded since their energy is of the order of V. 

Let us consider two consecutive sites and therefore, 
nearest neighbors of each other. Considering the sub- 
space with 1 = 0, there are three possible configurations 
for this pair of sites, which we will call links and they are 
{hp); (ph); (hh) where p stands for an occupied site and 
h for an empty one. The total number of these links in 
the AB2 chain is given by 

N,,p + Nph + Nhh = 4Afc (23) 

and Nfip = Nph- Note that unlike the case of the t-V 
chain, the number of links Nhp and Nhh iu the strong 
coupling AB2 chain is not a conserved quantity. 

A. Basis 

Let us consider a unit cell of the AB2 chain (which 
has three sites. A, B and C). The set of states for this 
cell correspond to five possible configurations in what 
concerns the particles distribution : i) zero particles with 
all sites being unoccupied; ii) one particle which may be 



at site A,B or C. The states with just one particle either 
at site B or site C can be combined to give a bonding 
and an anti-bonding state [see Fig. [l2|b) and Fig. [T2j[c)]. 
The anti-bonding state is a localized state as discussed in 
previous sections. The bonding state is a itinerant state; 
iii) two particles which must be at sites B and C in order 
to avoid a nearest neighbor occupied pair of sites [see 
Fig. [T2|a)]. These particles are also localized particles 
because if they hopped, a state with a nearest neighbor 
occupied pair of sites would be created. 

Let us consider now the case of N particles in a AB2 
chain with Nc unit cells. Let us define the operators 
Bl, = {Bj + Cj)/V2 and B^^ = {B\ - c\)/^. Note 
that the product of these two operators creates two par- 
ticles, one at site B^ and the other at site as expected. 
Using this new basis, the Hamiltonian can be rewritten 
in a simpler form 

PoifPo - -%/2t^e'^°/2^^ ( n X (24) 

+ n "^x4,J + l^.+2 +H.C. 





FIG. 12. (a) Two-particle localized state and one-particle 
(b) itinerant and (c) localized state involving B and C sites. 
The filled and empty circles indicate occupied and unoccupied 
sites respectively. Gray circles indicate bonding {+) and anti- 
bonding (-) occupations. The particles in (a) are localized due 
to the interaction {V/t — oo) and in (c) due to geometric frus- 
tration. In (b), the particle is free to hop to the neighboring 
A sites. 




FIG. 13. (a) State with localized as well as itinerant particles; 
(b) Open boundary region where the itinerant fermions in the 
state given in (a) may move. Note that the open boundary 
regions terminate always at B and C sites, (c) Corresponding 
state of the linear lattice. 

the following, we shortly review the main results of the 
algebraic solution. 

B. Itinerant states 



One could be tempted to say that the localized states play 
no role in this simplified Hamiltonian since only the -B^ ^ 
operator appears, but that would be incorrect. A basis 
for the ground state subspace {I = 0) can be constructed 
from states of the form 



Nbc 



Na 



I*) = n Bi^i n n n 4jo) 

i—l j — ^ n—1 m— 1 

(25) 

where the sets {a}, {13}, {7}, {^}, {a±l}, {/3±1}, {7±1} 
and {/i ± 1} have no common elements in order to sat- 
isfy the I — condition and Nbc, Nb_, Nb^ and Na, 
are respectively the number of localized pairs of fermions 
at sites B and C, of antibonding localized fermions, of 
itinerant bonding fermions and of itinerant fermions at 
A sites, so that 2Nbc + Nb_ + Nb^ +Na = N. Due to 
the projection condition, the localized fermions created 
by the j and ■B^_ - operators act creating open 
boundaries for the hoppings of fermions. The Hamilto- 
nian mixes states with different configuration for the A^- 

and b\_ j products but the part relative to the localized 
fermions remains always the same. In Fig. |13[ we illus- 
trate this. 

The diagonalization of the Hamiltonian given by 
Eq. p4|) is achieved by first distributing the localized 



fermions (if any) and then solving the remaining prob- 
lem for the itinerant fermions which move in regions 
confined by the localized particles, leading to a fac- 
torized form of the eigenstates, jitinerant fermions) iX) 
I localized fermions). The eigenvectors and eigenvalues 
of the ID t-V model for twisted boundary conditions as 
well as open-boundary conditions can be obtained by use 
of the Bethe ansatz, but an equivalent but simpler alge- 
braic solution is known in the strong-coupling limit. In 



Let us consider first the case when no particle is local- 
ized. In this case, the model becomes equivalent to the 
t-V chain apart from a renormalization of the hopping 
constant {t — )■ ^/2t). The solution of the t-V chain relies 
in recognizing that the condition I = leads to a conser- 
vation of the number of links N^p and N^h so that the 
tight-binding term only exchanges the position of these 
links^^. Interpreting the (hp) links as non-interacting 
particles hopping in a chain where the empty sites are 
the (hh) links, the solution is attained. Since the total 
number of {hh) and {hp) links \s L = L — N , the effec- 
tive chain is reduced in relation to the real t-V chain. 



Note that L is the number of sites in Fig. 13 c), after the 
mapping to a ID chain. The fact that the tight-binding 
particles occupy two sites of the real t-V chain leads to 
a twisted boundary condition which is dependent on the 
momenta of the tight-binding particles and the eigenval- 
ues are given by^ 



N 



E{{k}, P) = -2^21 ^ cos fc^ - - 



P 



L 



(26) 



with k = n-2'K /L, and P = n-2'n/L, with fi = Q, . . . , L—l 
and n = 0, . . . , L — 1. The set of pseudo-momenta {fc} 
and P must satisfy the following condition 



N 



P-^='^k (mod27r). 

^ ^=l 



(27) 



The mapping of this solution of the t-V chain (with even 
number of sites) into the AB2 chain without localized 
particles is direct, with odd sites corresponding to A sites 
and even sites to sites B and C (which will be unoccupied 
or in a bonding configuration). 
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Localized states 



Let us consider now a state where localized particles 
are present. Then one has one or more open boundary 
regions that terminate always in B and C sites as shown 
in Fig. 13 i.e., the number of sites is odd (the and 
Ci sites count as one site). The solution for the itinerant 
particles in one open-boundary region is rather simpl^^. 
Again, a mapping to a system of free fermions in a re- 
duced linear lattice is possible in a similar way, thinking 
of (hp) links hopping in a backg round of (hh) links, but 
it is simpler to state as irpsMl ^^j^g^^ ^j^g positions i of 
particles in the reduced chain are given by the relation 
i = i — Ni where i is the position of the particle in the t-V 
ID lattice and Ni is the total number of particles between 
the initial site and site i. The number of sites of the re- 
duced chain isi^L — A^+1. The energy contribution 
of these itinerant fermions is Ea = — 2-^2^ ^^j^ cos(fci) 
with k — 7rn/(iVrcd + 1) and n = 1,2, . . . , N^ed- Note 
that in the same state, and according to the distribution 
of the localized fermions, several confined regions may be 
present. 



D. Ground state 

The contribution of localized fermions to the energy of 
a given state is zero, but that does not necessarily im- 
ply that the ground state will always correspond to the 
minimum number of localized fermions. As an example, 
let us consider p — 1/3. An eigenstate for this filling is 
n^odd ^m|0)- The respective energy is zero since this 
state corresponds to having the itinerant fermion band 
completely filled, so that the positive kinetic energies bal- 
ance the negative kinetic energies. In this case, it is pos- 
sible to construct lower energy eigenstates with localized 
fermion pairs (created by the operator i3|c|), which al- 
low room for the itinerant fermions to move and therefore 
lower the total kinetic energy. 

This competition between itinerant and localized 
states will in fact start to occur at lower filling, p = 2/9, 
since at this filling, the reduced chain is half-filled and 
therefore positive kinetic energy states will start to be 
filled when adding additional itinerant particles. In or- 
der to lower the energy, one wants to keep the number of 
itinerant fermions just below half-filling and to have the 
maximum possible length for the reduced lattices where 
the itinerant particles move. Noting that a localized pair 
forbids the presence of particles in the two neighboring A 
sites, one concludes in order to maximize the number of 
sites available for the itinerant fermions these localized 
pairs should gather in a single cluster. 



The ground-state energy is obtained from Eq. ( 26 ) for 
fiUings less than 2/9. If A'^ is odd, all single-particle states 
with pseudo-momentum k between zL2tt/L - {N—1)/2 are 
occupied and ^ A: = 0. Therefore, 



odd 
GS 



= -2\/2t- 



™ (2JVe-Af) 



COS 



2N, 



(28) 



If N is even, all states with k between — 2Tr/ L-{N— 2) /2 
and 27t/L ■ N/2 or between -27r/Z • N/2 and 2^/1 ■ {N - 
2)/2 are occupied and ^ = ivr/L • N/L. So, 



sm 



= -2V2t- 



sm 



{ ttN \ 
\2N^-N J 

\2Na-N J 



TT ~ (f> 
2Nr_ 



(29) 



In the thermodynamic limit, the difference in the last two 
expressions becomes irrelevant and one can write 



N, 



2V2t 



(2 — 3p) sin 



np 



(30) 



For filling larger than 2/9, it is energetically favorable 
to have localized pairs of fermions in consecutive unit 
cells so that only one region exists for itinerant fermions, 
with 2Nc — 2Nbc ~ 1 sites and the respective reduced lat- 
tice will have L = 2Nc — N sites. The number of itinerant 
fermions is A^itin — N ~ 2Nbc- The number of BC local- 
ized pairs is such that the band for the itinerant fermions 
in the respective reduced lattice is as near as possible 
from half-filling, so that Nbc = Int(3A^/4 - NJ2) + 1 
or Nbc = Int(3A^/4 - A^c/2). In the thermodynamic 
limit, one can write A^itin = — N/2 and one has 
a closely packed localized cluster with length equal to 
Nbc/2 = 3N/8 - iV,/4. 

The ground-state energy in this case is 



tAoc 



-2V2t( 



TT 

2 2N, 



A^itin + 1 



iV+ 1 



sin 


7r 




2 




sin 


TT 


1 


2 





which for a large AB2 chain simplifies to 

2V2t^ 



A^. 



-(2-3p), 



(31) 



(32) 



where p = N/3Nc. Since this ground state is localized, 
no magnetic fiux dependence is present. Note that for 
A^ = 2A'c — 1 the reduced lattice has only one site and 
the only possible energy is zero. Also for A^ = 2Nc, one 
has a Wigner crystal like state with Nbc = Ac and zero 
energy. Also for A^ = 2Nc — 2, the reduced lattice has two 
sites and and the number of itinerant fermions is zero or 
two, therefore the only possible energy is zero. 

The density in the compact cluster of BC localized 
pairs present in the ground state for p > 2/9 has den- 
sity equal to 2/3, while the region available for the 
movement of the itinerant fermions has constant density 
A^itin/(|A^itin — 1) which for large A^itin is approximately 
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FIG. 14. (a) Ground state energy as a function of filling 
of the AB2 chain in the strong-coupling limit V — 00 and 
with Nc = 6. The curves are the analytical results given by 



Eqs. (281, (291 and 1 31 1 and the dots are the energy levels 
obtained by numerical diagonalization of the V = <x AB2 
chain, (b) Ground state energy in the thermodynamic limit 
as a function of filling in the strong-coupling limit = 00. 
The transition between a metallic and an insulating ground 
state occurs exactly at p = 2/9 and the minimum energy is 
obtained when p ~ 0.2. 



constant and equal to 2/9. So as one increases the filling 
in the AB2 chain, one is reducing the length of a phase 
with lower density and increasing the size of the higher 
density phase. 

Phase separat ion h as been observed in other ID and 
2D Hamiltonian^^SliU jj^ ^j^g j^q jnodel, phase sepa- 
ration occurs for J/t = 2.5 — 3.5, with the system divided 
into an electron-rich region and a hole-rich region^. An 
open-boundary Heisenberg model rules the dynamic of 
spin degrees of freedom of the electron-rich phase. Par- 
ticularly relevant to our study of the AB2 chain is the 
phase separation of the anisotropic Heisenberg model 
or XXZ model, which can be mapped onto the ID t- 
V model. In this case, phase separation corresponds to 
the appearance of the ferromagnetic phase in the XXZ 
model"^^. The anisotropy constant of the XXZ model, 
A = J\\/ J±_ under the Jordan- Wigner transformation be- 
comes A — —V/2\t\ and for b oth models the phase sep- 
aration occurs for A > llSillll^ that is, phase separation 



occurs for attractive nearest-neighbor interaction (ferro- 
magnetic interaction). In the t-V AB2 chain, we have 
shown that phase separation occurs for strong repulsive 
interaction. Note that a spinless fermion in the t-V model 
corresponds to an up-spin in XXZ model while a hole 
corresponds to a down-spin, and therefore the lower and 
higher density phases of the AB2 chain correspond to 
different magnetization regions of the antiferromagnetic 
XXZ AB2 model. 



E. Luttinger liquid description 

One of the most interesting points concerning the t-V 
AB2 chain is the following: is the t-V AB2 chain a Lut- 
tinger liquid? That is, can the low energy excitations of 
this model be described as bosonic charge den sity f luctu- 
ations governed by the harmonic HamiltoniarPlHll 

HLL^vsY.llHbg + ^ivNiN-Nof + vjA'! (33) 

In the LL Hamiltonian, N, J, L, vn, vj, and vs are, re- 
spectively, the particle number, current number, system 
length, particle velocity, current velocity, and sound-wave 
velocity. One is easily tempted to calculate these Lut- 
tinger liquid parameters from the previous strong cou- 
pling results, obtaining for large L and for p < 2/9, 



VN 



1 d^iEr/L) 36V2t 



TT d{N/LY (2 - 3p)= 



sm 



_d'^{ET/L) _ 9(2 -3p)t 



dWLy 



2V2 



9V2t 



sm 



2-3/9 

'Snp 
2-3p 



2- 3p 



inp 
2-3p 



(34) 



(35) 



(36) 



where P = x and vn, vj, and vs are respectively 
the particle, current, and sound velocities. Note that 
N/L = 3p/2 since L = 2A'c. These are results similar to 



those obtained for the strong coupling t-V ringSSEHin the 
thermodynamic limit. The Luttinger liquid parameter 



1 



vn/vs 



4 

(2 - 3p)2 



(37) 



determines the anomalous correlation exponents and is 
filling dependent reflecting the reduction of the effective 
size of the chain with filling. 

For filling larger than 2/9, the ground state has the 
localized pairs of fermions in a compact cluster so that 
only one region exists for itinerant fermions, with 2A'c — 
2Nbc — 1 sites and the respective reduced lattice will 
have L — 2Nc — N sites. One again is tempted to de- 
scribe the low energy behavior of this system of itinerant 
fermions in this open-boundary lattice using as above 
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a LL description. With open boundaries, one loses the 
translation invariance but it has been shown by several 
authors that an open-boundary bosonization can still be 
carried oulP^lEH. Many of these studies of LLs with open 
boundaries were motivated by the fact that the introduc- 
tion of a single local impurity in a Luttinger liquid breaks 
the ID system in two parts, at leas t in w hat concerns the 
low energy properties of the systenP^I^. There is a clear 
analogy between the role of these local impurities and 
the localized fermions in the AB2 chain. For very large 
L, one expects the solution of the system to become less 
dependent on the boundary conditions and therefore the 
LL velocities are the same as those of an equivalent sys- 
tem with periodic boundary conditions at the same filling 
and with the same lattice size (since the LL exponents are 
dependent on filling) . But we have seen that the filling in 
this open-boundary region is constant and equal to 2/9, 
and since the LL velocities depend only on filling, one 
obtains 



VN 



243t 



= 3V2t, 



vs 



27t 
27I' 



-2tp _ 



(38) 

Note that this ground state is N^. degenerate due to the 
translation invariance of the phase separation bound- 
aries. 

The missing point in the previous analysis is that 
the LL Hamiltonian does not describe the excitations of 
ground state that involve the creation of an additional 
localized anti-bonding fermion or an additional localized 
BC pair of fermions. However this is only important if 
these excitations are low energy states. In a qualitative 
picture, the strong coupling AB2 chain can be described 
as if there were a flat band at zero energy with Nc two- 
level sites. This is similar to the noninteracting case, but 
now the number of localized fermions in such flat band is 
twice as many. For p ^ 2/9, the creation of a localized 
anti-bonding fermion or a localized BC pair of fermions 
implies an excitation energy of the order of t and for en- 
ergies and temperatures much less than this value, the 
LL description is valid. As the filling approaches 2/9, 
these excitations become low lying (since the Fermi level 
goes to zero) and must clearly be taken into account in 
the low temperature description of the system. In this 
situation, even at very low temperature, the low energy 
AB2 set of eigenvalues becomes a complex mix of the sets 
of eigenvalues of LLs with different sizes, fillings, bound- 
ary conditions, and LL velocities. Despite these remarks, 
one should note that the strong coupling LL velocities do 
indeed characterize the correlations that do not involve 
the creation of additional anti-bonding states or addi- 
tional localized BC pairs of fermions. For example, a 
two-point ground-state correlation involving only A sites 
such as the Green's function restricted to A sites. 

For p > 2/9, the compressibility 



1 



1 d^Eoip) 



(39) 



K L dp'^ 

is infinite due to the linear behavior of the ground-state 



energy. This is the expected behavior of a phase sepa- 
rated ground state and is known to occur also in the t-J 
model as well as in the XXZ modeP^. In these mod- 
els, the compressibility diverges as (A^ — A)~^ as one 
approaches the phase separation critical value (which is 
Ac = 1 in the case of the attractive t-V chain) by in- 
creasing the interaction constant!^. The compressibility 
can also be calculated from the LL relation 



2K„ 



(40) 



In the XXZ model (or in the equivalent attractive t-V 
model), the compressibility diverges since the sound ve- 
locity vanishes and the Luttinger parameter Kp diverges 
as A — >■ Af,. From the last two relations, one easily con- 
cludes that as the filling in the AB2 approaches 2/9, the 
compressibility does not diverge, since the curvature of 
the plot of E{p) in Fig. [l4](b) does not vanish just be- 
low p = 2/9 or equivalently, neither the sound veloc- 
ity vanishes nor the Luttinger parameter Kp diverges as 
2/9. 

For p < 2/9, the curre nt vel ocity gives the charge stiff- 
ness at zero temperatur^^Eol 



2 a(0/i)2 



0c=O • 



(41) 



For p > 2/9, the charge stiffness is zero reflecting the 
open boundaries condition for the itinerant fermions and 
consequent zero dependence of the energy levels on mag- 
netic flux. This discontinuity of the Drude weight at the 
phase boundary is known to occur in other ID models 
with phase separation^^ . In this case, the current veloc- 
ity should correspond to a very small frequency peak in 
the optical conductivity since it is known that the open- 
boundary condition shifts the spectral weight associated 
to the Drude peak to a flnite frequency pealf^. 



VII. IMPLICATIONS FOR THE EXTENDED 
HUBBARD AB2 MODEL 

In this section, we discuss the relevance of the results 
obtained in this paper for the t-V AB2 model for the 
spinful extended Hubbard model in the AB2 geometry 
and in the strong-coupling limit U ^ t and U ^ V. 

First, let us recall the known facts about the strong 
coupling Hubbard model in a rin^^t^ and in a ID chain 
with open-boundary conditions^. The extended Hub- 
bard Hamiltonian for a ring with L sites is given by 

H = -i ^(c|^Ci+io- + 4+i^Ci^) 

i 

+U'^n.i^nii + V'^nini+i, (42) 

i i 

where the creation (annihilation) of an electron at site i 
with spin a is denoted by c]^ {cia) with Uia being the 
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number operator Uia = cl^Ci^ and n, = + ni^. When 
t — 0, all states with the same number of doubly oc- 
cupied sites and the same number J2aa' ^o-a' of nearest- 
neighbor occupied sites are degenerate, where N^^a' is the 
total number of nearest-neighbor pairs with spin config- 
uration aa' . The eigenvalues of the extended Hubbard 
model in the atomic limit are given by E{Nd, {Na-a'}) = 
■ U + J2cra' ^o-o-' • V. Here we will only address the 
low energy subspace of the strong-coupling limit where 
Nd = to make this discussion simpler. 

If we consider the Hubbard ring with t <C U, but 
V ^ 0, one has the so-called Harris-Lange model^ and 
the model eigenfunctions can be written as a tensorial 
product of the eigenfunctions of a tight-binding model of 
independent spinless fermions in the ring (where the spin- 
less fermions are the electrons deprived of spin) and the 
eigenfunctions of an Heisenberg model (with exchange 
constant J = t^JU) for the spins of the electrons in a re- 
duced chairPSmSl j'Yie spinless fermions ring is threaded 
by a fictitious magnetic flux <f> — generated by the spin 
configurations in the reduced Heisenberg chain (where qs 
is the total spin momentum) and the eigenvalues to order 
t are given by 



E{ki,. . . , fcAT^+AfJ = 2t ^ COS (jti 



i=l 



(43) 



where kt = {2tt /L)ni, rii = 0, . . . , L — 1 are the momenta 
of the holes in the spinless ring. 

The nearest-neighbor interaction is obviously indepen- 
dent of spin of the electrons that occupy nearest-neighbor 
sites and it is easily introduced in the previous picture 
so that ioT U ^ t and U ^ V, the Nd = spectrum of 
the Hubbard model to the order of t is that of interacting 
spinless fermions in a ring threaded by a fictitious mag- 
netic flux 4> = qs (since the nearest-neighbor repulsion 
between electrons leads to an nearest-neighbor repulsion 
between spinless fermions), that is, the extended Hub- 
bard model in a ring in the limit [/ —> oo has the energy 
dispersion of the spinless t-V chain. The only effect of 
spin in this limit is the generation of a fictitious flux and 
the increase of degeneracy. 

If now we consider the strong coupling Hubbard chain 
with open-boundary conditions, the solution is even sim- 
pler. As stated in Ref.l^, for p < 1, the physics of the 
model is the same as that of a spinless tight-binding 
model. In fact, the electrons hop along the chain but 
the spin configuration of the electrons remains always 
the same and becomes irrelevant for the determination 
of the eigenvalues of the model. The different spin con- 
figurations only contribute to the degeneracy of the en- 
ergy levels. The same reasoning as above can be followed 
and one concludes that the extended Hubbard chain with 
open-boundary conditions and in the limit U ^ t and 
U ^ V has the same energy spectrum as the t-V chain 
with open-boundary conditions. 

Let us now consider the extended Hubbard model in 



the AB2 geometry, 

H = Ho + vY,{^f + nf+i) {nf + nf) 
i 

+ U J2{nf^nf^ + nf^nf^ + ngn^), (44) 

i 

where 

AT, 

Ho = -tJ2Il [e''"^''^^ i4<^B,, + bIA,+i.) (45) 



a j = l 



H.c. 



We consider the limit J7 — > 00 where no doubly occupied 
sites are possible. 

In the AB2 geometry, we can not state as we did for 
the Hubbard chain that for y = the model can be 
mapped into a system of spinless fermions in a tight- 
binding AB2 chain threaded by fictitious magnetic fiux 
generated by the spin configuration momentum. The rea- 
son is that in the case of the Hubbard chain, only hop- 
pings at the boundaries generate a different spin configu- 
ration and all the different configurations generated can 
be obtained from one of them applying a circular per- 
mutation operatoi^^^Hll!, However, in the AB2 geometry, 
besides the circular permutation of the spin configuration 
due the hoppings at the boundaries, additional spin mix- 
ing can also be generated by electrons hopping along one 
plaquette (which circularly permute a subset of two or 
three spins). However, this additional exchange process 
requires the existence of an occupied nearest-neighbor 
pair of sites in that plaquette as an intermediate step. 
Therefore, as V is increased, this exchange process is in- 
hibited. 

Let us therefore address the extended Hubbard AB2 
model in the limit t <^ V <^ U. In this case the addi- 
tional plaquette exchange does not occur. We show below 
that the low energy spectrum of this model, that is, the 
set of eigenvalues corresponding to eigenstates with no 
doubly occupied sites and no occupied nearest-neighbor 
pair of sites will be the same as that of the t-V AB2 
model in the limit t ^ V, but with enlarged degeneracy 
just like for the extended Hubbard ring in the strong- 
coupling limit. The picture behind this result will be a 
mix of the pictures presented for the Hubbard ring and 
for the Hubbard chain with open-boundary conditions. 

First we show that the itinerant and localized states 
of the extended Hubbard AB2 model remain the same 
as those of the strong coupling t-V AB2 model in what 
concerns the charge distribution, but now one has a con- 
siderable larger set of states due the spin degrees of free- 
dom. We adopt the same basis as that of the previous 
section, but now spin must be considered, that is, the 
bonding and anti-bonding states have a spin degree of 
freedom. In the previous section we have shown that two 
type of localized fermions could occur: i) one-particle lo- 
calized states corresponding to a single particle in one 
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plaquette in a standing wave state, which obviously re- 
mains locahzed if the particle has spin; ii) two-particle 
BC localized states induced by the large nearest-neighbor 
repulsion. The spinless localized pair maps into four lo- 
calized BC pairs with different spin configurations which 
again remain localized due to the large nearest-neighbor 
repulsion. Concerning the itinerant particles, the Hamil- 
tonian does not alter the spin of the bonding states as 
the particles hop along the AB2 chain. 

Therefore, we can construct the eigenstates of the ex- 
tended Hubbard AB2 model in the limit case t <^ V <^ U 
in an equivalent way to that followed in the previous sec- 
tion. If no localized particles are present, the extended 
Hubbard AB2 model in the limit case t <^ V U can 
be mapped into the extended Hubbard ring in the same 
limit and the energy levels are given by 

N / 

E{{~k},P)^-2V2tJ2cosfk-^-^^] (46) 

i—l^ ^ 

where hi and L have the same definition as for the t-V 
ring in the strong-coupling limit, has the same defini- 
tion as for the Hubbard chain and is a external mag- 
netic flux (but keeping a zero flux in each plaquette). In 
the thermodynamic limit, the effect of the fictitious and 
external fluxes becomes irrelevant. 

If localized particles are present, they create open 
boundary regions where itinerant electrons (in bonding 
states when at sites B and C) will move according to 
the extended Hubbard AB2 model in the strong-coupling 
limit. Each of these open boundary systems has a fixed 
number of particles and can be mapped as above into the 
extended Hubbard chain in the same limit, but now with 
open-boundary conditions. In this stated before, 

the spin configuration is irrelevant and contributes only 
to the degeneracy of the energy levels which are exactly 
those described in the previous section for the spinless t- 
V AB2 model with open-boundary conditions. Since the 
energy levels in the thermodynamic limit are the same as 
those of the t-V AB2 model in the strong-coupling limit, 
phase separation will also occur in the extended Hubbard 
AB2 ring and precisely at the same filling. This exact cor- 
respondence required the strong-coupling limit and was 
obtained only for the low energy subspace. Higher en- 
ergy subspaces involve the presence of occupied nearest- 
neighbor pairs of sites or of double occupancies, but a 
similar approach can in principle be followed again with 
a mapping to higher energy subspaces of the extended 
Hubbard chain. 

For intermediate values of the nearest-neighbor inter- 
action y, the mapping is no longer exact, but qualita- 
tively the features described for the t-V AB2 chain should 
be expected in the extended Hubbard AB2 chain (with 
[/ 3> i and [/ ^ y). For example we expect a simi- 
lar evolution for the ground-state energy with decreasing 
V , but with different evolution of filling. Note that the 
ground-state energy for non-interacting electrons (with 
spin) as a function of filling has the same form of Fig. [5] 



but the spin degeneracy of the single-particle eigenvalues 
implies that the filling values as well as the ground-state 
energy must be multiplied by two. This is equivalent 
to distributing spinless particles among two independent 
AB2 chains. 



VIII. CONCLUSION 

As mentioned in the introduction, the ID anisotropic 
Heisenberg model (the XXZ model) can be mapped using 
the Jordan- Wigner transformation into the ID spinless t- 

V model (with an additional "chemical potential" term, 
which can also be interpreted as an on-site energy). The 
ground state filling of the t-V model is related to the 
ground state magnetization of the XXZ model, which 
can be controlled by the application of an external mag- 
netic field to the XXZ model. This external field does 
not change the eigenstates of the XXZ model, but cre- 
ates an additional chemical potential term in the t-V 
model which allows to control the filling of the ground 
state. A ferromagnetic Heisenberg interaction means 
that the spinless fermions attract each other while an an- 
tiferromagnetic exchange implies a repulsive interaction 
between nearest-neighbor spinless fermions. A Jordan- 
Wigner transformation is also possible in quasi-lD or 2D 
lattices (see Ref.'^ for a review), and in particular, the 
XXZ AB2 model in the strongly anisotropic limit can be 
mapped into the strong coupling t-V AB2 model with 
additional phase factors which are nonlocal (and there- 
fore usually treated in a mean-field approaches. In this 
paper, we have not gone into the details of the transfor- 
mation of the XXZ AB2 model into the t-V AB2 model 
(which would require the determination of the phase fac- 
tors and the introduction of the on-site energy terms), 
but some conclusions can be drawn on general argu- 
ments. First, one-magnon localized states in the strongly 
anisotropic XXZ AB„ model corresponding to standing 
waves in a t-V AB„ array will exist since no additional 
phase factors appear when only one spinless fermion is 
present (only one spin flip). Also, states with several 
spin flips corresponding to localized magnons in differ- 
ent plaquettes are also eigenstates since the XX term 
of the XXZ model gives a zero contribution in the re- 
gions between the localized magnons (since all spins are 
aligned). So the one-particle localized states in the t- 

V AB2 model correspond to localized and independent 
magnons created in a ferromagnetic background, which 
have been observed in frustrated magnetic systems un- 
der high magnetic fields but below the saturation fielcP^. 
Second, independently of the phase factors, the two- 
particle localized BC particles will remain the same be- 
cause the nearest-neighbor interaction obtained in the 
Jordan- Wigner transformation retains the same form as 
in the ID case, and since these two particles are com- 
pletely localized, the phase factors are irrelevant. Fur- 
thermore, these localized particles create open bound- 
ary regions for the itinerant fermions independently of 
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the phase factors. So, we expect a similar behavior of 
the XXZ AB2 model in the strongly anisotropic limit 
to the one described in this paper for the strong cou- 
pling t-V AB2 model. Such two-particle localized states 
in the strong coupling t-V AB2 model correspond to lo- 
calized pairs of magnons in the the XXZ AB2 model in 
the strongly anisotropic limit. 

The results presented in Sec. IV for single-particle 
states in AB„ lattices may also be relevant to Joseph- 
son junction arrays (with the same geometry) both in 
the quantum limit and in the high capacitance (classi- 
cal) limit. It is known that a Josephson junction AB2 
chain with half a flux quantum per plaquette exhibits 
a highly degenerate classical ground state reflecting the 
completely flat energy bands of the AB2 tight-binding 
model for this value of flujP^. The tunneling of Cooper 
pairs between the different superconducting islands in 
an AB2 geometry can also be described using a bosonic 
tight-binding AB2 model. Furthermore, in the high- 
capacitance limit, the charging energy due to the elec- 
trostatic interaction within each (large) superconducting 
island can be neglected and the Hamiltonian becomes a 



classical XY model, H = —JJ2{ij) cos 
where is the superconducting phase of the island i, J 
is the Josephson coupling and Aij is the phase shift due 
to the presence of an external magnetic field, obtained 
from the integral of the vector potential along the path 
from i to j. This Hamiltonian can be mapped into a one- 
particle tight-binding model with the same geometry and 
under the same magnetic flux^^. In most geometries, the 
minimum energy phase configuration of the Josephson 
junction array will be obtained from the state of mini- 
mum energy of the tight-binding model (if this state is 
homogeneous). In the most general case, the Hamilto- 
nian can be interpreted as the mean energy of a phase 
vector and the stable phase configuration of the Joseph- 
son junctions array is obtained minimizing this energy 
and this minimization may imply mixing tight-binding 
states of different bands. 

Let us now compare the exact results of the strong cou- 
pling AB2 chain with the mean-field results. Comparing 
Figs. 11 and 14 (b), one concludes that the mean-field 
results overestimate the interaction energy contribution 
for p < 2/3. Basically, within the mean-field approach 
the low energy band increases its energy with increas- 
ing V, missing the fact that itinerant states are pos- 
sible, which avoids the positive energy contribution of 
the nearest-neighbor interaction as found in the strong- 
coupling limit. The mean-field results are qualitatively 
correct in what concerns the energy interval for the 
ground-state energy for p < 2/3. In Fig. 11 



one can 



see that the ground-state energy remains negative even 



for large V in the density range < p < 2/3. The min- 
imum of the ground-state energy as a function of filling 
shifts continuously to lower filling in contrast with the 
V = 00 result where such minimum occurs at p ~ 0.2. 
Also the large slope for p > 2/3 agrees with the fact that 
for such fillings nearest-neighbor pairs are present and 
such states have infinite energy when V = 00. 

To conclude, in this manuscript, the spinless AB2 
chain with nearest-neighbor Coulomb interactions has 
been studied for any filling and taking into account mag- 
netic flux. In the case of independent fermions, a sim- 
ple construction of the localized states that generate the 
flat bands both in the presence and absence of flux has 
been found and generalized for ID or 2D arrays of quan- 
tum rings. The V/t versus filling phase diagram of the 
AB2 chain was obtained using a mean-field approach. 
The dependence on filling of the mean-field ground-state 
energy agrees qualitatively with the exact ground-state 
energy for infinite V. The ground-state energy for in- 
finite nearest-neighbor repulsion has a quantum criti- 
cal point at filling 2/9 where a metal-insulator transi- 
tion occurs. This transition refiects the phase separa- 
tion between a high density phase {p = 2/3) and a low- 
density phase {p — 2/9) that occurs at fillings larger 
than 2/9. Such phase separation occurs because the in- 
finite nearest-neighbor repulsion leads to the appearance 
of two-particle localized states (besides the one-particle 
localized states due to the topology of the AB2 chain). 
These localized states create open boundary regions for 
itinerant carriers and in order for these itinerant fermions 
to have only negative kinetic energy, phase separation be- 
comes favorable. At low filling, the low energy properties 
of t-V AB2 chain can be described by the spinless Lut- 
tinger Hamiltonian, but for filling near or larger than 
2/9, the AB2 set of eigenvalues becomes a complex mix 
of the sets of eigenvalues of LLs with different fillings, 
boundary conditions, and LL velocities. If the itinerant 
fermions have spin, but a very strong on-site repulsion 
is present (that is, in the case of the extended Hubbard 
model in the strong-coupling limit C7 ^ 1/ ^ i), the 
energy-dispersion relation to the order of the hopping in- 
tegral remains the same as that of the spinless AB2 model 
in the presence of a flux and phase separation occurs at 
the same filling. 
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